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Abstract: When analysing dark matter halos forming in cosmological n-body simulations it is common 
practice to obtain the density profile utilizing spherical shells. However, it is also known that the systems 
under investigation are far from spherical symmetry and rather follow a triaxial mass distribution. In 
this study we present an estimator for the error introduced by spherically averaging an elliptical mass 
distribution. We systematically investigate the differences arising when using a triaxial density profile 
under the assumption of spherical symmetry. We show that the variance in the density can be as 
large as 50% in the outer parts of dark matter halos for extreme (but still credible) axis ratios of 
0.55 : 0.67 : 1. The inner parts are less affected but still show a scatter at the 16% level for these 
prolate systems. For more moderate ellipticities, i.e. axis ratios of 0.73 : 0.87 : 1, the error is smaller 
but still as large as 10-20% depending on distance. We further provide a simple formula that allows to 
estimate this variance as a function of radius for arbitrary axis ratios. We conclude that highly prolate 
and/or oblate systems are better fit by analytical profiles that take into account the triaxial nature of 
cosmological objects. 

Keywords: cosmology: theory — cosmology: dark matter — methods: analytical — methods: n-body 
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1 Introduction 

The current standard model of structure formation, in 
which the Universe is dominated by cold dark matter, 
predicts that the radial density profile of dark matter 
ha los is best descr i bed b y the functional form proposed 
bv lNavarro et alJ (Il997ft (hereafter NFW): 
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with the two fitting parameters p c and r s , and a was 
found to be equal to unity, i.e. a = 1. 

During the last decade great strides have been un- 
dertaken to verify and better understand the origin 
of this density profile of dark matter halos forming in 
cosmological n-body simulations. The functional form 
proposed by NFW is independent of cosmology and 
the two fitting parameters p c and r s are actually re- 
lated. Equation Q therefore describes a one-param eter 
family of functions. While INavarro et alJ fl997) still 
relied on systems containing of the order 10 4 particles 
(and hence resolving the profile down to about 10% of 
th e virial radiu s, acco rding to the convergence criteria 
of iPower et alJ (l200.?t) ') studies containing about two 
orders of magnitude more particles within the virial 
radius were only able to confirm their findings 



Fukushige fc Makinoll997l:lFukushige fc Makinol200„ 
Moore et al.ll998l:lJinal2000UJing fc Sutol2000HJing fc Sut 

2002). rlowever, there has been a debate about the ex- 
act value of the log arithmic inner slope ranging from 
a ~ 1 to 1.5 (e.g.. INavarro et al.l[l997l. iMoore et alJ 
1999). Recently a great deal of effort has gone into 
even higher resolution simulations which have revealed 
density profiles of cold dark matter halos down to scales 



120051 : IPiemand. Moore fc Stadell2004l) . The highest 
resolved simulation to-date reached an effective mass 
resolution of about 130 million particles in a cluster 
sized dark matter halo still supporting the evidence 
for a central c usp with a logarithm ic inner slope of 
about a ~ 1.2 (Diemand ct al. 2005). However, these 
steep central density cusps are not supported by obser- 
vations. By measuring the rotation curve of a galaxy, 
one can in principle infer the density profile of its dark 
matter halo. High resolution observations of low sur- 
face brightness galaxies are though best fit by halos 
with a core of constant density (Simon et al. 2005; 
Swaters et al. 2003; de Block & Bosma 2002). This 
discrepancy between the prediction of theory and ac- 
tual observations has caught the attention of cosmol- 
ogists for nearly a decade now and is far from being 
resolved. 

Aside from addressing the mass distribution of dark 
matter halos from a simulator's point of view, a fair 
amount of labour went into a theoretical understand- 
ing of the origin of the NFW profile equation But 
only when making strong assumptions about a connec- 
tion between the density and the velocity dispersion 
can one can make analytical predictions. For example, 
under the assumption that the phas e-space density is a 
power law in radius as suggested bv lTavlor fc Navarro! 
12001T) one can solve the spherical Jeans equation to 
Tibtain the density profile iPehnen fc McLaughlinl2005l . 
lAustin et all200 5. Hans enl2004h . These studies, guided 
by the results of the numerical simulations, confirm the 
inner logarit hmic slope to be i n the range a ~ 1 to 2. 
But see also Hansen fc Stadell ll2004t) for recent claims 
that equilibrated dark matter structures have a central 
0.8. 



well below one percent of the virial radius ||Fu kushige. Kawai^^Maicino 

l2004l : lTasitsiomi et all2004UNavarro et alj2004ilReed^t^n It needs to mentioned that (nearly) all these stud 
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ies are based upon the assumption of spher ical symme- 
try ut ilizing radially distributed shells. But lJing fc Sutol 
( 2002) already noted that halo density profiles as found 
in cosmological simulations are better fit by apply- 
ing an elliptical or triaxial modeling of dark matter 
halos. Detailed investigations of the distributions of 
shapes of halos forming in cosmological n-body sim- 
ulations now demonstrate that the majority of ob- 
jects has intrinsic shapes best characterized by axis 
ratios of about 1:0.74:0.64 fe.g-lKasun fc Evrarcl2005l: 



the surface of a sphere of a given radius due to the ellip- 
ticity of the halo. There are hence parts on the sphere 
with values larger than the spherically averaged den- 
sity as well as parts with lower densities. However, we 
though refer to a 2 (r) as the "error" measure (or error 
estimate) as it quantifies the mean deviation (either 
positive or negative) from the actual (triaxial) density 
distribution under the assumption of spherical symme- 
try. 



iBailin fc Steinmetzl200Rl . lAllgood et all2006ft . lOeuri et ail 

(2005), for instance, showed that such prolate mass 
distributions can have a significant influence on the in- 
terpretation of lensing signals observed in galaxy clus- 
ters, i.e. halo triaxiality may act to bias the profile 
constraints derived under the assumption of a spher- 
ically symmetric mass distribution. It therefore ap- 
pears mandatory to systematically study the error in- 
troduced by the supposition that halo profiles are spher- 
ical while their actual distribution is triaxial. Can the 
debate about the central cusp be connected to this er- 
ror? Answering this question is the major motiviation 
for the current study. 




2 The Equations 

A model for a triaxial NFW den sity profile is given by, 
for instance, Uing fc Sutol (120021) and reads as follows: 



g(R) 



(2) 



R/R a (l + R/Rs) 2 ' 
where the relation between the spherical and elliptical 
coordinates for a given Cartesian point x,y, z 



r/R s 



Figure 1: The error estimate equation (thick 
lines) for the most likely axis ratio a : b : c (solid 
line) and the 1 — a scatter (dot-dashed and dashed 
lines) as derived from cosmological n-body simula- 
tions. The faint horizontal lines correspond to the 
error when using an isothermal density profile as 
given by equation JJJ). 



x = r sm a cos < 
y = r sin 9 sin <j 
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(3) 



depends on the axes a < b < c of the ellipsoid and 
hence the degree of triaxiality T (e.g. iFranx et alJ 
1991) 



(4) 



If we now average the triaxial density g(R) on the 
surface of a sphere with radius r 



(Q(r)) r 



1 
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g(R)sm9d6d<f> , (5) 



where g(R) = g(R(r,8, <f))) is a function of r, and <j> 
through equation we can obtain a measure for the 
error for the variation of the actual density g(R) by 
calculating the dispersion 




Figure 2: Variation of fitting parameter A with 
a/c along lines of constant triaxiality T. 
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in 



g(R) - (g(r)) 



(e(r)) 



sin 6d6d<f) 



(6) 

We need to stress that this quantity is per defini- 
tion the typical dispersion in the density for points on 



3 The Results 

In figure Q we now plot the error estimate as given by 
equation © (thick lines) for a number of axis ratios. 
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Figure 3: Same as figure El but showing (3 as a 
function of a/c. 



where A and /? are free parameters. The best fit curves 
to the error equation © are also shown in figure as 
thin solid lines. 

The question now rises about the dependency of 
the fitting parameters A, /3 on the axis ratios a : b : c 
and the triaxiality T, respectively. To this extent we 
show in figures [5] and [3] how A and /3 vary as a func- 
tion of a/c along lines of constant triaxiality T. We 
note that the variation in both parameters is rather 
small and the average value of the median values for 
constant T amount to (.4) = 0.017 and </3) = 0.530. 
When adopting these parameters, equation JH) pro- 
vides a (universal) formula to estimate the error in- 
troduced by spherically averaging an elliptical mass 
distribution for given axes ratios a/c and b/c, respec- 
tively. 

4 Conclusions 



The ratios chosen agree with the find ings of cosmologi- 
cal w-b ody simulations as reported bv lKasun fc Evrardl 
(2005). Their distributions are close to Gaussian and 
we chose to present equation JSJ for axis ratios based 
upon the mean ratios (a/c) and {b/c) as well as the la 
scatter about these values. This leaves us with the ra- 
tios (a/c) = 0.64±0.09 and (b/c) = 0.77±0.10. Hence 
we plot equation @ for the ratios 



• a : b 

• a : b 

• a : b 



c = 0.64 : 0.77 : 1 (mean), 
c= 0.73 : 0.87 : 1 (+la), and 
c= 0.55 : 0.67 : 1 (-la). 



While the errors are rather small in the central 
parts we note a discrepancy of up to 50% in the outer 
regions of dark matter halos for the more triaxial axis 
ratios. 

To better understand and gauge the functional be- 
haviour of the curves presented in figure Q we calcu- 
lated the error estimate equation @ for the same axis 
ratios but using a simple isothermal density profile 



g(R) = 



Qc 



(R/Rs 



(7) 



The results can be viewed in figure Q too, and are 
represented by faint horizontal lines for the respective 
axis ratios. We note that for a single power-law density 
profile the error is constant with radius. We further 
observe that the crossing point between the isothermal 
and the NFW profiles occurs about the scale radius 
(i.e. r/R a ~ 1) at which the NFW profile approxi- 
mates the isothermal behaviour. This suggests to fit 
the radial dependence of the NFW error estimate by a 
curve with a functional form similar to the logarithmic 
slope of the density profile itself, i.e. d log g/d\og R. 
We found the following function to accurately describe 
the data 



Motivated by the findings of cosmological simulations 
that dark matter halos have triaxial mass distributions 
Ce.g-.lKasun fc Evrardll2005l IrBailin fc Steinmetzl2005l . 
Alleo od et, alJl2006l) we investigated the errors in the 
derivation of the density profile under the assumption 
of spherical symmetry. This supposition is rather com- 
mon and nearly all results presenting evidence for the 
"cuspy" behaviour of dark matter halos rely upon this 
premise. While others already showed that a spherical 
assumption will lead to discrepanc ies in the interpre- 
tation of observa tional data (e.g., iJing fc Sutoll2002l : 
lOguri et alJl2005Tl we presented an estimator to quan- 
tify this error. 

We set out to shed light on the question whether 
the discrepancy and debate abou t the exact value of 
the logarithmic cen tral slope fe.g. jNavarro et alll997l . 
iMoore et al.lll99 9f) can be related to the error intro- 
duced by the simplifying assumption of spherical sym- 
metry. We now conclude that the inner parts are only 
marginally affected and the error will be too small to 
explain variations in the slope as large as ±0.5. A 
change in a from, for instance, 1.5 down to 1 entails 
a difference in the central density of a factor of about 
3 (at approx. 10% of the scale radius) rather than 
the low percentage level indicated by equation © and 
seen in figure Q However, there are claims that dark 
matter halos tend to be more elliptical in the central 
regions (iDubinski fc Carlbergll99l1.|Jing fc Sutol200 ' 
| Kazantzidisl2004lBailin et alkoOSUAllgood et all2Qi 
iPower et al.ll200a) . But those shape changes are only 
minor and will hence not be able to lift the error to 
the required level. 

We further provided a fitting function that allows 
to estimate the error introduced by spherically averag- 
ing a triaxial mass distribution as a function of radius 
for (almost) arbitrary axis ratios a/c and b/c. 
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